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Abstract. Physical universality of a cellular automaton was dehned 
by Janzing in 2010 as the ability to implement an arbitrary transfor¬ 
mation of spatial patterns. In 2014, Schaeffer gave a construction of a 
two-dimensional physically universal cellular automaton. We construct a 
one-dimensional version of the automaton. 
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1 Introduction 

A cellular automaton (CA) is a finite or infinite lattice of deterministic finite 
state machines with identical interaction rules, which, at discrete time steps, 
update their states simultaneously based on those of their neighbors. They are an 
idealized model of massively parallel computation. From another point of view, 
the local updates can be seen as particle interactions, and the CA is then a kind 
of physical law, or dynamics, governing the universe of all state configurations. 

We study the notion of physical universality of cellular automata, introduced 
by Janzing in [1], which combines the two viewpoints in a novel and interesting 
way. Intuitively, a cellular automaton is physically universal if, given a finite 
subset D of the lattice and a function h on the shape-1? patterns over the states 
of the CA, one can build a ‘machine’ in the universe of the CA that, under the 
dynamics of the CA, decomposes any given pattern P and replaces it by h{P). 

A crucial point in this definition is that we need to perform arbitrary compu¬ 
tation on all patterns, not only carefully constructed ones. This has quite serious 
implications: The machine M that takes apart an arbitrary pattern of shape D 
and replaces it by its image under an arbitrary function h is not in any way 
special, and in particular we are not allowed to have separate ‘machine states’ 
and ‘data states’ with the former operating on the latter. Instead, we can also 
think of M as a pattern, and must be able to construct a larger machine M' 
that takes M apart and reassembles it in an arbitrary way. 
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This notion differs essentially from most existing notions of universality for 
CA such as intrinsic universality [7] , universality in terms of traces as discussed 
by Delvenne et al. in [3], and the more well-known I-know-it-when-I-see-it type 
of computational universality promoted by Wolfram in [5]. In these notions, 
one can usually implement the computations and simulations in a well-behaved 
subset of configurations. Physical universality bears more resemblance to the 
universal constructor machines of Von Neumann [B] , which construct copies of 
themselves under the dynamics of a particular cellular automaton, and were the 
initial motivation for the definition of CA. Another property of CA with a similar 
flavor is universal pattern generation as discussed in [2, meaning the property 
of generating all finite patterns from a given simple initial configuration. 

In Janzing’s work [4] some results were already proved about physically uni¬ 
versal CA, but it was left open whether such an object actually exists. A two- 
dimensional physically universal CA was constructed by Schaeffer in [5], but it 
was left open whether this CA can be made one-dimensional. We construct such 
a CA, solving the question in the positive. 

2 Definitions 

Let A be a finite set, called the state set, and a grid. A cellular automaton 
is a map / : A* —>• A* defined by a local function F : —>• A, where 

N C is a finite neighborhood, so that /(x)i, = F^xn+v) for all x G A* 
and D G It is reversible if there is another CA g : A* —>■ A* such that 
fog = gof = id. An example of a reversible CA is the shift by n gIA, defined 
by cr"'(x)„ = Other examples can be constructed as follows. Let A be a 

Cartesian product n^=i with projection maps tt^ : A —>• A^, let ni,..., G 

be arbitrary, and let 7 : A —>• A be a bijection. Then the CA / defined 
by /(x)i, = ^{TTi^x-v+ni), ■ ■ ■ j'^kixv+rik)) is reversible. We call / a partitioned 
CA, and the components Ai are tracks. In the CA, the tracks are first shifted 
individually by the rii, and then the bijection 7 is applied to every cell. 

A CA / is physically universal if the following condition holds. For all finite 
domains D,E C Z'^, and all functions h : A^ A^, there exists a partial 

configuration x G A^ and a time t G N such that for all P G A^, we have 
f{xUP)E = h{P). It is effectively physically universal if t is polynomial in the 
diameter oi D U E and the computational complexity of h according to some 
‘reasonable’ complexity measure. In this article, we use circuit complexity, or 
more precisely, the number of binary NAND gates needed to implement h. One 
could reasonably require also that the configuration x is computed efficiently 
from the efficient presentation of the function h. Our proof gives a polynomial 
time algorithm for this. See Section [ 8 ] for a discussion. 

3 The Cellular Automaton 


Our physically universal automaton is a partitioned CA / defined as follows. 


— The alphabet is A = {0,1}^. 

— The shifts are 2, 1, —1 and —2, and we denote S = {2,1, —1, —2}. 

— For each a,b € {0,1} bijection 7 maps the state (1, a, b, 1) to (1, b, a, 1), and 
(o, 1,1, b) to ( 6 ,1,1, a). Everything else is mapped to itself. 

Intuitively, in the CA / there are four kinds of particles: fast right, slow right, 
slow left and fast left. At most one particle of each kind can be present in a cell. 
On each step, every slow (fast) particle moves one step (two steps, respectively) 
in its direction. After that, if two fast or two slow particles are in the same cell, 
then the direction of every particle of the other speed is reversed. This resembles 
the two-dimensional CA of Schaeffer, where particles move to four directions 
(NE, NW, SE, SW) with speed one, and the head-on collision of two particles 
causes other particles in the same cell to make a u-turn. 
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Fig. 1: A sample spacetime diagram of /, and a schematic version on the right. 
Time increases upward. Particles are represented by arrowed bullets. For exam¬ 
ple, a bullet with arrows to the east, northeast and northwest represents three 
particles moving at speeds 2 , 1 and — 1 , respectively. 


4 The Logical Cellular Automaton 

For the proof of physical universality, we define another CA on an infinite alpha¬ 
bet. For this, define the ternary conditional operator by p(a, b, c) = c©(aA(c© 6 )) 
for all a,b,c G {0,1}. That is, p{a, b, c) is equal to 6 if a = 1, and to c otherwise. 
In many programming languages, p{a, &, c) is denoted by a ? & : c. 

Definition 1. Let V = be an infinite set of variables, and denote 

by T the set of Boolean functions over finitely many variables ofV. The logical 
extension of / is the CA-like function / : -G A^ on the infinite alphabet 

A = , where the four tracks are first shifted as in f, and then the function 

(a, b, c, d) I—>■ {p{h A c, d, a),p{a A d, c, b),p(a A d, b, c),p{b A c, a, d)) 


is applied to each coordinate. A valuation is a function v :V ^ {0,1}. It extends 
to T and then into a function v : —>■ A^ in the natural way. 






The logical extension simulates multiple spacetime diagrams of /: one can 
see that the definition of / is equal to that of /, except that each particle is 
replaced by a Boolean formula that corresponds to the conditional presence or 
absence of a particle. We think of >1 as a subset of A containing the constant 
0 or constant 1 function in each track. Note that / is also reversible, and we 
denote by f~^ its inverse function. See Figured] for a spacetime diagram of /. 
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Fig. 2: A schematic spacetime diagram of /, where a, /? G Particles are repre¬ 
sented by solid lines, and Boolean particles by dotted lines. Note that Boolean 
particles can be created, even though / itself conserves the number of particles. 


The following result holds basically by construction. 

Lemma 1. Let x G he a configuration, and let v : V ^ {0,1} be a valuation, 
so that v{x) G AA'. Then for all t G Z we have f*{v{x)) = v{f*{x)). 

The idea of the proof of physical universality of / using this new CA is the 
following. We may assume that D — E = [0,n — 1] in the definition of physical 
universality, for some n G N. Then, we construct a spacetime diagram of / with 
the following properties. First, in the initial configuration x G A^, the cells of the 
interval [0, n — 1] contain 4n distinct variables from V. All other cells of x contain 
either 0 or 1. There also exists t > 0 such that /*(^)[o,n-i] contains the Boolean 
functions computing the function h in the definition of physical universality. In 
the course of the construction, we define which cells of x contain a 1. 

Definition 2 . We introduce the following terminology for the construction. 

— The configuration of interest, denoted by x G initially contains the fully 
general’ state ( 04 ^, 04 ^+ 1 , 04 ^+ 2 , a 4 i+ 3 ) in every cell i G [0,n— 1], and 0 
everywhere else. During the construction, we change the cells of x from 0 to 
1, but keep referring to it as x, so some of the definitions below depend on 
the stage of the proof. The spacetime diagram of interest is defined similarly. 



— A (spacetime) position is an element of 1 a x S (1 x 1 a x S, respectively), 
representing a (spacetime) point that may contain a particle of certain speed. 
Note that time is bi-infinite, since our cellular automata are reversible. 

— There is a Boolean particle at spacetime position {i,t,s) if TTs(j(x)i) is not 
the constant 0 function, and a particle if it is the constant 1 function. 

— There is a collision at position {i,t) if f*{x)i contains at least three Boolean 
particles, and a crossing if there are at least two Boolean particles. 

— The input is the pattern X[o^„_i] € A". 

— The gadget is the contents of x outside [0,n — 1], an element o/. 

— A line is a subset of 1 a xN of the form L — L{i, t,s) = {(* + st, t) | < € N} 
for some speed s € S and i € Z. It is occupied (in a time interval I C Z) if 
one of its coordinates (in the region Z x I) contains a crossing or a Boolean 
particle of speed s. We denote = i st and L* = {L^^\t, s). The set of 
occupied lines in the spacetime diagram of interest is denoted Cocc- 

For example, there are three crossings in Figure [21 two of which are collisions. 
The highest intersection of two dashed lines is not a crossing, as it does not take 
place at an actual coordinate (one of the white circles). Every line segment in 
the figure defines an infinite occupied line. 


5 The Diffusion Lemma 

As stated above, we initialize the gadget to the all-0 partial configuration, in 
which situation we have the following lemma. It states that any finite set of 
particles in the CA / eventually stop interacting and scatter to infinity. The cor¬ 
responding result is proved for the physically universal CA of [5] by considering 
an abstract CA over the alphabet {0, i, 1}, with the interpretation that ^ can 
be either 0 or 1, and this lack of information is suitably propagated in collisions. 
In our version, the role the new state ^ is played by Boolean particles. 

Lemma 2 (Diffusion Lemma). Let x G be such that xt = 0 for all i ^ 
[0, n—1]. Then there are 0{n^) crossings in the two-directional spacetime diagram 
of x under f that all happen in a time window of length 0{n). For all other times 
t G Z, there are 0{n) Boolean particles in f^{x), and they are contained in the 
interval [—2|t|, n -|- 2|t|]. 

Proof. We prove the claim in the positive direction of time. By induction, one 
sees that after any t > 0 steps, there can be no right-going Boolean particles in 
the cells (—oo, t — 1], and no left-going particles in the cells [n — t, oo). After these 
sets intersect at time |"n/2], there are no collisions, so no new Boolean particles 
are created. Thus the number of Boolean particles is at most 6n, that is, twice the 
length of the segment of (x) that may contain Boolean particles. We may 

also have 0{n^) crossings between Boolean particles going in the same direction 
with different speeds. Thus there are 0{n^) crossings in total. □ 



Fig. 3: An illustration of Lemma [2J The dashed line is the configuration x, and 
the thick segment is the interval [0,n — 1]. All collisions take place in the dark 
gray area, and Boolean particles may occur in the light gray area. The eight 
horizontal line segments are the scattering Boolean particles, grouped by speed. 


At this point, in the configuration of interest, we have an empty gadget, 
and the spacetime diagram contains 0{n) Boolean particles at any given time. 
Since the CA / is reversible, the values of the corresponding Boolean expressions 
determine the values of the original variables. 

6 Manipulating the Spacetime Diagram 

6.1 Controlled Modifications 

In this section, we introduce new particles in the gadget that will collide with 
the existing Boolean particles and create new ones. This is called scheduling 
collisions. We never schedule a collision on an occupied line, and never add a 
Boolean particle on an existing crossing. This is formalized in the following. 

Definition 3. A modification of the gadget is (a, 6 , C, T, t, /)-controlled for num¬ 
bers a, 6 G N, sets of lines C and T, time t £ 1, and interval I G 1, if the 
following conditions hold: 

1. the modification consists of adding at most a particles to the gadget, 

2. at most a new occupied lines and b new crossings are introduced, 

3. no existing crossings become collisions, 

4- no line in LUCocc is occupied by a new Boolean particle or a new crossing, 

5. no spacetime position in the set 

F{T, t, I) = {(^, t, s')\i£l,s' £D}\r (1) 

gets a new Boolean particle, and 

6. no line in T gets a collision after time t. 














If the conditions hold in the time interval (—cx),i] (in particular, condition 6 
need not hold at all), then the modification is weakly (£, T, /)-controlled. 

In practice, a controlled modification is one where we add to the gadget a 
hnite number of particles that affect the spacetime diagram of interest only where 
we want it to be affected: the spacetime positions in T*. The lines in CUCocc and 
the positions near T, that is, those in F{T, t, I), are ‘protected’ from accidentally 
obtaining any auxiliary Boolean particles created in the modification. 

Definition 4. Let j,t € Z. The positive cone rooted at {j,t) is the set of space- 
time coordinates C{j,f) = {(j +/,t+ t') | t' € N,j' € [—t',t']}. 

The following lemmas are parametrized by the numbers mi for i = 1, 2, 3,4. 
They denote, intuitively, the number of existing crossings, the number of occu¬ 
pied lines, the size of the additional ‘forbidden area,’ and the number of particles 
involved, respectively. Also, the expression ‘if t' = l7(iV(mi, m 2 , m 3 , m 4 )) is large 
enough, then P{t') holds’ for t' € N, a function > N and a property 

P means that there exists a number T < K ■ /(mi,m 2 ,m 3 ,m 4 ) -I- K', where 
K,K' > 0 are constants independent of the m^, such that P{t') holds for all 
t' > T. 

6.2 Moving the Boolean Particles 

We now prove that we can add a finite number of particles to the gadget, so that 
after some number of steps, a collection of Boolean particles is ‘moved’ onto any 
desired lines, at the same time. We do this one particle at a time, and without 
interfering with the trajectories of any other existing particles. 

Lemma 3. Suppose we have a spacetime diagram of interest with mi crossings 
and m 2 occupied lines, and a Boolean particle p with label (3 € P on a line 
Lp that contains no collisions after time t. Let C be a collection of m 3 lines 
not containing Lp. Let L ^ C be an unoccupied line that passes through some 
spacetime coordinate {j',t') € C(j,t) with t' > t. Let T be a set of 0{ms,) lines 
containing L, and let I Ch be an interval of length Oljn^). Ift' = t-\- f2{mi + 
m 2 -I- m 3 ) is large enough, then there is an {0{l),0{m2),C,T,t',1)-controlled 
modification after which the spacetime position L* contains a Boolean particle 
with label (3. The same holds if the line L is unoccupied only in the time interval 
{—oo,t'], but the modification is weakly controlled. 

Proof. Denote the speed of the particle p by s, and the target speed by s'. We 
assume for simplicity that s = 2 and s' = —2. One can perform a similar analysis 
for the cases (s, s') = (2, —1) and (s, s') = (1, —2). By mirror symmetry, and by 
repeating such a modification at most twice, one obtains all possible movements. 

We change the Boolean particle p to a slow right Boolean particle after 
some k gN steps, by scheduling a suitable collision at the spacetime coordinate 
(/ -\-2k,t-\-k). After some k' gN more steps, we turn this Boolean particle into 
a fast left one; this happens at the spacetime coordinate (j 3-2k 3-k',tk 3-k'). 




Our goal is that t' — t — k — k' steps after the second collision, the Boolean particle 
is at Of course, whether this happens depends on our choices of k and 

fc'. Namely, we must have 

= j + 2k + k' — 2(t' — t — k — k') = j + 4:k + 3k' + 2t — 2t' 
due to the speeds of the Boolean particles. This implies that 

fc'= i(j'- j + 2(t'-t)-4fc) (2) 

so that in particular f — j + 2{t' — t) — Ak should be divisible by 3, or in other 
words, k = Ak = j' — j + 2{t' — t) mod 3. 

The other obvious constraints are 0 < k, 0 < k' and k + k' — Expand¬ 
ing ([ 2 ) in the latter two inequalities, we obtain 

max(0, / -j - {t' -t))<k< -j + 2(t' - t)) (3) 

for the variable k. Denote by / C Z this interval of possible values for k. 

By the assumptions, we have j — (t' — t) < j' < j + (t' — t), which implies 
f — j — {f — t) <0, so that the left hand side of ([3]) can be replaced by 0. The 
length of the interval I is therefore 

a= - {j — j' 2 {t' — t)) > -{t' — t) = Q{mi +m 2 + m 3 ). 

Since the only other requirement on k is parity modulo 3, there are 17 (toi +m 2 + 
TO 3 ) possible choices for k. If we have k G I and k = j' — j + 2(t' — t) mod 3, we 
say k respects the spacetime constraints. 


Now, for any k respecting the spacetime constraints, we give a set of particles 
that, when added to the gadget at time 0, realize these two collisions assuming 
that they do not collide with any existing particles. The positions (coordinate- 
speed tuples) of the particles are Pi = {j — 2t — k', 2), P 2 = {j + 3k + t, —1), 
P 3 = (j +4:k + 2t, —2) and P 4 = {j + t + 3k + 2k', —1). We denote by Li the line 
containing the position Pi at time 0, and also Lp = Lq. 

First, let us forget about the lines C and the existing Boolean particles other 
than p, and track how p (continued backward to its imaginary starting position 
on X, although it may of course have been produced in a collision), and the new 
particles added at Pi, P2, P3, P4, interact. We show that, in fact, they move the 
particle p to as desired for any choice of k that respects the spacetime 

constraints, and all in all, no new collisions other than the two desired ones 
are introduced. Of course, once we have done this, it is enough to show that 
for some choice of k that respects the spacetime constraints, the conditions of 
Definition [ 3 ] are satisfied, since then the new Boolean particles introduced by 
our modification do not interfere with the existing ones in unwanted ways. 

Let us proceed to the analysis. The reader should consult Figure 0 ] for a 
visualization. At the spacetime coordinate {j -I- 2 k, t -f k), the Boolean particle 
p collides with the particles added to P2 and P3, resulting in all four types 
of Boolean particles. The new occupied line we are interested in is the path 
where we are moving p, or the slope -1 line L5. At the spacetime coordinate 
(j + 2 k + k',t + k + k'), where P3, P4 and the slowed-down p collide, all four 
Boolean particles are produced as well. Again, the interesting line is the one 
where we move p, or the line L = Lq with slope — i. 

Let us compute the equations of the lines Lq, Li,..., Lq: 

Lo = {{j -2t + 2y,y) \ y G N}, Li = {{j - 2t - k' + 2y,y) \ y G N}, 

L 2 = {{j + t + 3k - y,y) \ y G N}, L 3 = {{j + 2t + Ak-2y,y) \ y G N}, 

Li = {{j + t + 3k + 2k' -y,y) \y gN}, L 3 = {{j - t + k + y,y) \ y G N}, 

Le = {(/ + 2 t' - 2 y,y)\y G N}. 

Now, let us compute their intersection points (in addition to the desired inter¬ 
sections (j -I- 2 k, t + k) and (j + 2 k + k' ,t + k + k')) into Table [T] The last two 
intersections, D and E, are not shown in Figured since D did not fit, and E 
does not appear in the configuration (the others always appear). 

It is easy to verify from the equations of the lines that for each i G { 1 , 2 , 3 , 4 , 5 }, 
every choice of k results in a different line Li (the lines Lq and Lq are always 
the same). Also, the two collisions we introduced are different for each k. Recall 
now that there are mi crossings and m2 occupied lines in the spacetime diagram 
of interest, and that the set C has cardinality m3. Furthermore, the cardinal¬ 
ity of the set T of target lines and of the set E{T,t',I) of forbidden spacetime 
positions defined in (IT|) is likewise 0(1713). Thus, there are 0 (mi + m2 -\- m3) 
such choices for k that at least one of the lines Li is in £ U £occ or contains an 
existing crossing, a line in £ U £occ gets a new collision, or a spacetime position 
in F(T, t, I) gets a Boolean particle. If a = l 7 (mi -I- m2 -I- m3) is large enough. 


Table 1: The incidental intersections of the lines Lq, Li,..., Lq. 


Lines 

Intersection point 

Label 

Pi and P 3 

[j -1- 2 fc - \k',t-\- k Ifc') 

A 

Pi and P 2 

(j-t 2 fc- ifc') 

B 

pQ and P 4 

{j + 2 k + ^k',t k ^k^) 

C 

P 2 and Pg 

{j "1“ 2/i: + ^k\t k ^k') 

D 

P 3 and P 4 

[j -|- 2 fc -|- ik', t k — 2k') 

E 


we can guarantee the existence of a A: that respects the spacetime constraints, 
and is not in the aforementioned set of cardinality 0 {mi + m 2 + m 3 ). 

For such k, the modification of adding the particles Pi, P 2 , P 3 and P 4 to 
the gadget satisfies the conditions 3, 4 and 5 of Definition [3] in the time interval 
(— 00 , t']. Condition 2 is also satisfied in this interval, since the only new occupied 
lines are Li,..., Pg) and every new crossing happens at the intersection point of 
an existing occupied line and one of Pi,..., Pg, the number of which is 0 ( 7712 ), 
or at one of the five intersection points in Table [TJ Finally, condition 1 is satis¬ 
fied by definition, and thus the modification is weakly (0(1), 0 ( 7772 ), P, T, t', /)- 
controlled, and results in the spacetime position (j', t', —2) containing a Boolean 
particle with label /3. 

If the target line P is completely unoccupied, then condition 3 is satisfied 
in the entire spacetime diagram, so that no new collisions or occupied lines are 
introduced after time t'. We can also check from Tabled] that the only crossing 
on line P is at {j + 2 k + k\ t + k + k'), so that condition 6 is satisfied. Then the 
modification is (0(1), 0 ( 7772 ), P, T, t' , /)-controlled. □ 

The point of the set of forbidden lines P and positions I is that we can in fact 
move an arbitrary number of Boolean particles at the same time. The idea is 
that we move them by applying Lemma |3| repeatedly to one Boolean particle at 
a time, always adding the target lines of the remaining ones into the protected 
set P. This guarantees that the lines are not accidentally occupied too early. 

Corollary 1 . Suppose we have a spaeetime diagram of interest with mi eross- 
ings and m2 occupied lines, and some Boolean particles pi, ■ ■ ■ ,Pmi with labels 
Pk & ^ on lines Lp,^ that contain no collisions after time t. Let C be a eolleetion 
of m3 lines not containing any oftheLp,^. Let Lk ^ C be unoccupied and mutually 
disjoint lines that pass through some spacetime coordinates {jk,t') G C{Lpl,t) 
with t' > t. Denote D = {Pi, ■ •., Pm4}, ond let L G Z be an interval of length 
0(7713). If t' = t + 0(7771 + W3 -|- 7774(7772 + 7714)) is large enough, then there 
is an (0(7774), 0(7774(7772 -\-m4)), C,T,t', I)-controlled modification after which 
each spacetime position L\, contains a Boolean particle with label Pk- The same 
holds if the lines Lk are unoccupied only in the time interval (—oo,P], but the 
modification is weakly controlled. 











Proof. We apply Lemma[ 3 ]to each Boolean particle pk separately. First, let Ci = 
£ U {L2, ■ ■ ■, Lm^}. We can do a (weakly) ( 0 ( 1 ), 0 (m 2 ), £1, T, t,/)-controlled 
modification to the gadget such that the spacetime position L{ contains a 
Boolean particle with label / 3 i. This modification does not introduce crossings or 
Boolean particles to the lines of £i (in particular, it does not affect the Boolean 
particlesP2j • ■ • jPrm), and it creates 0(1712) new crossings and a constant number 
of new occupied lines. 

We continue inductively for k < m4, assuming that the particles pi, ... ,pk 
have been moved to their places by introducing 0(km2 + k'^) new crossings and 
0 (k) new occupied lines. We now define £k+i = £ U {Lk+2, ■ ■ ■ ,Lm}, and use 
Lemma| 3 ]to do a (weakly) {0(l),0(m2+k), £k+i,T, t, /)-controlled modification 
to the gadget, which causes to contain a Boolean particle with label f 3 k+i- 
Note that £k+i has size 0(1713 + 1714), there are 0 (mi + km2 + fc^) crossings 
and 0(7712 + k) occupied lines in total, and t' = t + Q(mi + m3 + 7714(7712 + 
7774)). Furthermore, if the lines do not intersect at all, no new crossings are 
introduced on them at any point. Thus, the conditions of Lemma [ 3 ] hold. 

After the inductive construction, we have modified the gadget by adding 
0(7714) particles to it. The total numbers of new crossings and occupied lines 
are 0(7714(7712 +7774)) and 0(7774), respectively. The remaining conditions for the 
modification to be (weakly) (0(7774), 0(7774(7772 + 7774)), £, T, t', /)-controlled are 
easy to verify. This concludes the proof. □ 


6.3 Computing with the Boolean Particles 

Next, we will do some computation with the Boolean particles. Namely, we show 
that the NAND of two Boolean particles can be computed nondestructively in 
any spacetime position, as long as we have enough time, and the target line is 
in the intersection of the cones rooted at the input particles. 

Lemma 4 . Suppose we have a spacetime diagram of interest with mi crossings 
and m2 occupied lines, and two speed-1 Boolean particles pi,p2 labeled jdi, (32 € IF 
with P2 £ Pi ^ 1 which occupy distinct lines Li and L2 that contain no collisions 
after time t. Let £ be a set of m3 lines not containing Li and L2, and let 
L ^ £ he an unoccupied line of slope 1 to the left of Li and L2 that passes 
through some spacetime coordinate (j',t') G C(Lf\t) D C(L2\t) with t' > t. 
Let T = {Li, L2, L}. If t' = t -\- L 2 (mi + m2 + m3) is large enough, then there 
is an (0(1), O (m2), £,T~,t',$)-controlled modification after which the spacetime 
positions L\ , L\ and L* contain Boolean particles labeled Pi, P2 and -<(PiAP 2), 
respectively. 

The condition /32 A/3i ^ 1 is there mainly for terminological reasons, ensuring 
that the output of the computation is actually a Boolean particle, and does not 
restrict the usefulness of the lemma, since the constant 0 function is available for 
us anyway. Similarly to the fact that any number of particles can be moved, we 
can now compute an arbitrary Boolean function, given enough time and space. 


Proof. We proceed as in the proof of Lemma[31 but skipping most of the details. 
For i £ {1, 2}, we will schedule a collision at the spacetime coordinate {ji+ki,t+ 
ki) for some ki < k 2 £ N, in which a slow left-moving and a fast right-moving 
particle collide with pi, creating a fast left-moving Boolean particle with label 
Pi and leaving pi intact. At time t + k 2 + k for some fc G N, we schedule a third 
collision that slows down the copy of pi. After some time, the still fast copy of 
P 2 will catch it, and we schedule a fourth collision at this point that creates a 
fast right-moving Boolean particle with label Pi A ^ 2 - A fifth collision will be 
scheduled at the spacetime position where this particle reaches the target line 
L, producing a slow right-moving Boolean particle with label ~'{Pi A P 2 ). 



Fig. 5: Computing the NAND of two Boolean particles. The auxiliary Boolean 
particles are shown in gray for clarity. 


Figure [S] shows the auxiliary particles needed for the collisions. It also shows 
that, for suitable values of /ci, k 2 and fc, the only collisions resulting from the 
modification are the five desired ones. As in Lemma[3l there are then a bounded 
number of new occupied lines and 0 ( 1712 ) new crossings, and if t' is large enough, 
we can find such values for the parameters that the conditions of Definition [3] 
are satisfied. □ 

Corollary 2. Suppose we have a spaeetime diagram of interest with mi colli¬ 
sions and m 2 occupied lines, a coordinate j £ Z, and Boolean partieles pi,... ,Pmi 
with labels Pk £ if at the spacetime positions (j + k — l,t,l). Suppose further 
that the lines Lp^, that the pk occupy do not contain collisions after time t. Let 
H : {O,!}’”^ —>■ {O,!}*” with m < 7714 be a Boolean function realizable with C 
NAND-gates such that o H ^0 for all £ = 1,..., m. Let C be a set of m 3 
lines not containing the Lp,^, let Li,... ,Lc+m ^ C be unoccupied lines of slope 
1 that pass through the spacetime coordinates (j — 1 ,t)(j — C — m,t), and 


let t' > t. Let also T = {Lc, ■ ■ •, Lc+m}- If 

t' = t + f2((C + TO 4 )(C'^ + mi + TO 3 + (C + m4)(m2 + m 4 ))) 

is large enough, then there is an ( 0 (C+m 4 ), 0{{C+m 2 +m 4 ){C+ 1714 )), C, T, t', 0)- 
controlled modification after which the spacetime positions contain Boolean 

particles with labels H(Pi, ..., (3rn)i- 

Proof. This proof is similar to that of Corollary[T] The C auxiliary lines are used 
for the computation of H one NAND-gate at a time. Let to = t, m^°^ = 0 and 
m® = 0 . 

For each i = l,...,C' + m, we apply Lemma 0] to some pair of particles pyi 
and pi ^2 on some of the lines Lp^ for fc e {1 ,..., m 4 } or Lk for k < i, and the 
target line Li, denoting the set of these three lines by Ti- We denote the number 
of new crossings after round z — 1 of the calculation by m^* and that of new 
occupied lines by We define Ci = £U{Li+i,..., L^+c}, so that the latter 

lines will not accidentally become occupied, and then \Li | = Olrn^) + C + m — i. 
For large enough 


ti = ti-i + n{mi + + m 2 + m 2 + m 3 + C + m — z), (4) 

there is an (0(1), 0 (m 2 T, 0)-controlled modification of the gad¬ 

get that causes the NAND of the particles and pp 2 to appear on the line Li 
at time U. 

When all the particles are in their places, it remains to bound the number of 
new particles, occupied lines and crossings in the spacetime diagram of interest, 
and the final time t' = tc+m- First, the number of particles and occupied lines 
increases by a constant each round, so we have m^*^ = 0(i). The number of 
crossings increases by m^*^ — m^* = 0 (m 2 + m^ = 0 (m 2 + i) on round 

i, since each new occupied line crosses each already occupied line at most once, 
so we have m^*^ = 0{m2i + z^). In particular, the total number of new crossings 
is = 0{{C + m 2 + m 4 ){C + m 4 )). Any large enough time difference 

ti — ti-i = l7(mi -I- m 2 Z + i'^ + m 3 -I- C -I- m 4 ) suffices (where we have plugged 
the values of m\ and m^ into o and simplified a bit), so we can choose 

tc+m - to = + m 4 )(C'^ -I- mi -I- m 3 -I- (C -I- m 4 )(m 2 -f m 4 ))), 

which finishes the proof. □ 

7 Physical Universality 

We can now prove our main result, the effective physical universality of /. 
Theorem 1. The cellular automaton f is effectively physically universal. 


Proof. We prove the result by a construction that uses the results from the 
previous sections. The reader should use Figure H] as a guide. Without loss of 
generality, let D = E = [0,n — 1], and let h : A" —)• be any function. We 

can view /i as a Boolean function H : {0, —>■ {0, on arrangements 

of particles in the domain D. We set the initial configuration of interest to 
be X S with the fully general pattern containing the 4n Boolean particles 
with labels ai,... ,a 4 „ S V in the input, and 0 in the gadget. Our goal is to 
modify the gadget so that for some time tgnai > 0 polynomial in n and the 
circuit complexity Ch of H, we have /*'‘”“-‘(x)i = H{ai,..., o: 4 n)i 4 i, 4 i+ 3 ] for all 
i € [0,n — 1]. Lemma [T] then gives the claim. 

In the construction, we apply Lemma [2] twice, first to the input pattern in 
the forward direction, to transform it to a dispersed pattern with 0{n) Boolean 
particles labeled Pi,, Pm € T, and then to the output pattern in the re¬ 
verse direction, to find a dispersed pattern with Oin) Boolean particles labeled 
(5i,..., bm’ G at positions (ji, Si),..., (jm', Sm') € Zi X S that / transforms 
into the output. The function H' : {0,1}™ —>■ {0,1}^ that we actually compute 
maps Pi,..., Pm to those elements of ^i,..., Sm' that are not identically zero, 
denoted , ■ ■ ■, Sn^., which is possible since the logical CA / is reversible. The 
function H' has circuit complexity C < Ch + 0{n'^), since it can be computed by 
simulating in reverse the OlrP) collisions that produce Pi,..., Pm from the in¬ 
put, then the function H, and finally in reverse the 0{n?) collisions that produce 
the output from Sn ^, •. ■, Suf. ■ 

The construction happens in five stages, scheduled at times 0 < tdis < tcoii < 
tcomp < tass < tfinai with large enough differences. As stated above, we will 
produce the dispersed pattern S ^,. •., Sn,, at time tass, and tfinai = tass + 0 {n) 
is the time when this pattern has evolved into the desired output. Let I = 
[—2(tfinai — ts,ss),'n'/4: + 2(tf}nai — tass)] be the interval containing the dispersed 
output pattern, and let be the size-O(n) set of all lines that pass through a 
coordinate of I at time tass- 

Diffusion Stage. In the Diffusion Stage, we simply apply Lemma [5] to the initial 
configuration of interest. It states that for some time tdis = 0{n), there are no 
crossings or collisions in the time interval [tdis,oo), and that the initial space- 
time diagram of interest contains 0{n) occupied lines and 0(n^) crossings. At 
time tdis, we have the aforementioned pattern of 0{n) Boolean particles, labeled 
Pi,..., Pn' € T, moving away from the origin at different speeds. 

Collection Stage. We now redirect the dispersed Boolean particles into a slow 
right-moving formation. For this, let tcoii = tdis + 0{C -I- n), and for i G Z, 
denote by Li = L(i,tco\i,l.) the line passing through position (i, 1) at time 
tcoii- Let £coii = {L-i I i = 1,..., C -I- fc} be the C + k lines reserved for 
computation, and let T^oii = {Li \ i = 0 ,...,m — 1} be the target lines for 
the inputs of the computation. Since there are 0{n) occupied lines and O(n^) 
crossings in the spacetime diagram of interest. Corollary [T] shows that there is 
an (0(n), O(n^), £coii T^oii, Aoii, 0)-controlled modification that moves the 

Boolean particles onto the lines Li,..., Lm at a large enough time tcoii- 


Computation Stage. In this stage, we compute the labels , • ■ •, (^n*. by sim¬ 
ulating the circuit representation of H'. For that, let tcomp = ^coii + 6*((C + 
n)^) and T^omp = {L-c-i \ i = l,...,fc}. Corollary [5] shows that there is an 
{0{C -\- n),0{{C + n)^), T^omp, ^comp, 0)-controlled modification which re¬ 
sults in k Boolean particles with labels Sm, ■ ■ ■, 5^, to be placed on the lines 
L-C-i, ■ • ■, L_c-k- 

Assembly Stage. In the Assembly Stage, we assemble the dispersed output 
pattern in such a way that no auxiliary particle affects its evolution before 
the final time tfinai- For that, let tass = tcomp + 6 *((C' -|- n)^), and let 7Iss = 
{L{jn.,tassT Sm) I i = l,...,fc} be the target lines that pass through the po¬ 
sitions of the Boolean particles in the dispersed output pattern at time tass- 
From Lemma [5] we know that the lines in TIss do not intersect in the time 
interval (—oo,tass]- We can then apply Corollary [T] again to obtain a weakly 
(0(n), 0(n(C + n)), 0, %ss, ^ass, /)-controlled modification that places a Boolean 
particle with label on each spacetime position {jrnTtass, Sm) while keeping 
the other spacetime positions in / x {tass} x S empty. 

Reverse Diffusion Stage. Here we no longer need to do any modifications. Since 
we explicitly forbade any auxiliary Boolean particles from entering the spacetime 
segment I x {tass}, the dispersed output pattern on it will evolve in tfinai — tass = 
0{n) steps into the final output pattern encoding H(ai ,..., a 4 „). □ 

8 Final remarks 

Our proof of the physical universality of / can readily be turned into a polynomial 
time algorithm that, given a circuit computing the function h : > A^ 

in the definition of physical universality, computes the corresponding gadget 
and the polynomially bounded number tfinai- However, this algorithm will need 
polynomial space, as it compares the new positions of auxiliary particles to all 
existing ones. In fact, for technical reasons that could be easily avoided, namely 
due to our choice of handling constant -0 particles separately, constructing the 
gadget is P-complete. To construct the gadget in logarithmic space, it might be 
necessary to fix particular choices of where the auxiliary particles are put. We 
have chosen the more abstract route in the hope that our methods generalize 
more directly to a larger class of CA. 

The existence of a physically universal CA was asked in [3] without fixing the 
number of states. Our CA has 16 states and radius 2. It would be interesting to 
find the minimal number of states and the minimal radii for physically universal 
CA. Of course, one can make any physically universal CA have radius I by 
passing to a blocking presentation, but this increases the number of states. From 
our CA, one obtains a physically universal radius-I CA with 256 states. 

Question 1. Are there physically universal CA on the binary alphabet? Which 
alphabet-radius pairs allow physical universality? 

A long list of open questions about physical universality is also given in [ 8 ]. 
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Fig. 6: A schematic diagram of the proof of Theorem[Tl not drawn to scale. 
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